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The baryonic decays of J/ip provide a new way to study the internal structure of 
' baryons. We apply a simple diquark model to the calculation of the decay cross-sections 

for the reactions J/ip — ► pp, p7V*(1440), N* N* , AA and E°E°. The results are different 
from those given by the ordinary constituent quark model. Hence these reactions may 
provide a new check of two different pictures for the baryons. 

^ . 

PACS: f4.20.Gk; f3.25.Gv 

Key words: decays of J/tj>; hadron production; diquarks 



<3\ 

o 
m 
o 



43 
Oh 



43 



- 1—1 

X 



C3 



1 Introduction 



In spite of a long history of baryon spectroscopy there are still many questions 
without clear answers. We know that the baryons are composed of three valence 
quarks, sea quarks and gluons but we are not certain if the valence quarks have 
constituent or current nature, if they cluster into diquarks or are well separated 
from each other, whether the sea quarks are in the form of qq soup or mesons etc. 
For a long time we have also missed direct sources of information on the properties 
of nucleon excitation states and could build our knowledge about them almost 
entirely on results from partial wave analysis of irN scattering data. The situation 
is changing dramatically with new experimental data coming from facilities such as 
CEBAF at JLAB, ELSA in Bonn, GRAAL in Grenoble or from BEPC in Beijing. 

A long-standing problem in N* physics is about the nature of the Roper res- 
onance AT* (1440) which is considered to be the first radial excitation state of the 
nucleon. However, various quark models have difficulty to explain its mass and elec- 
tromagnetic coupling, so it was suggested that it may be a gluonic excitation of 
the nucleon, a hybrid baryon. An ideal tool for studying the properties of N* reso- 
nances are the decays of J/ip — ► ppir and J/ip — > ppit+n^ since in these processes 
the tt°p and p7r+7r~ systems arc limited to isospin 1/2 states. The J/ip decays into 
baryon-antibaryon states discussed in the present paper provide another way to 
probe the internal structure of baryons. 

The J/ip decay cross sections for pp, pN* N*N* , AA and E°E° final states were 
calculated in Ref.p]|2] by using a simple quark model. In the present paper we use 
the diquark model and look into the possibility of forming the final state baryons 
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from the qq (generated in J/ip decay) and diquark-antidiquark (DD) pairs (cre- 
ated as vacuum excitation) . The model resembles the standard quark-pair-creation 
model S L used extensively to describe the mechanism of hadronic decays. Our di- 
quark model calculation gives a different prediction for these processes which can 
be compared with the previous calculation by the quark model ^12]. Future exper- 
imental data on these reactions may provide a check on these two different pictures 
for baryons. 

2 Diquark model formalism 

We assume that the decay of J/ ip meson leads to a creation of qq pair and the 
final state baryons are formed by their coupling to the diquark-antidiquark pair 
produced with the vacuum quantum numbers J PC = ++ . Restricting ourselves to 
the SU(3) flavour sector we have the following diquark nomenclature 0] 



(qq') r 


*> (ud — 


du) 


scalar-isoscalar diquark 


(qs) r 


(us — 


su) 


scalar-isodublet diquark 




uu 




vector-isotriplet diquark 


[qs] r 


- (us + 


su) 


vector-isodublet diquark 


[ss] r 


*> ss 




vector-isosinglet diquark, 



in which q (or q') stands for either the up or down quark and we use the round and 
box brackets for the axial scalar and the axial vector diquarks, respectively. The 
effective diquark-antidiquark vacuum creation amplitude can be written as 

(DD | V eS | 0) = f DD (-lf°- M °(T D M B T D M D \00) r, . (1) 

A similar structure was adopted in Ref. [5] for the DD — > qq process. The coupling 
parameter /q d may be different for the various DD species listed above. However, 
for the sake of simplicity we shall assume that this parameter does not depend on 
the diquark type. The phase r) = 1 stands for the creation of scalar diquarks and 
r\ = (— l) 1_cs <5 Cfl _ CD corresponds to the scalar coupling of vector diquark (with 
polarization ed) and vector antidiquark (with polarization ejj). 

The intermediate qq pair produced in J/tp decay can be either in the 3 5i or in 
the 3 Di state. In the J/ip rest system and for the initial J z = J — 1 polarization 
the relative qq wave function can be written as 

(p I qq, J = U z = 1) = clmMp)Ylm(P) | qq, (1 - M)) (2) 

LM 

where | qq, fi) stands for the spin- flavour part of the qq wave functions and p denotes 
the relative qq momentum. The radial wave functions (/>l(p) and the coefficients clm 
are normalized as 

El iT d P <fo(p) = 1 . Em c Im = 1 

(coo = 1, c 20 = v/T/lO, on = -v/3/16, c 22 = v/3/5). 
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Table 1. The spin-flavour overlap amplitudes [see Eq.@] derived for various BB final 
states. All coefficients V q (B) and the summed amplitudes f q V q are presented in units 

of f DD /2. 



BB state 


V U (B) 


V d (B) 


Vs(B) 


12f«v q 


pp,pN* ,N* N* 
AA 
E°E° 


1 - 1/(9^3) 


4/(9^2) 


-2/(9^3) 


4/(9 V2) 




2/3 
-2/(9^3) 


[1 - l/(3>/3)]/„ 
2/./3 

8/„/(9>/5) - 2/./(9v^) 



The pair couples with the D-D pair forming the final state baryon B and an- 
tibaryon £?. The spin-flavour overlap amplitudes follow from the SU(3) quark- 
diquark decomposition of the baryonic states (see Appendix). It is easy to show 
that 

(BB, fx' | V cS | qq, fx) = V q (B) , (3) 

where we introduced the BB spin triplet states in the same fashion as the qq 
spin-flavour states. The amplitudes © depend on both the baryon specification 
as well as on the flavour of the intermediate quark. For a reference we give the 
amplitudes relevant for the considered processes in Table [3 It is worth noting that 
the amplitudes J3J vanish if they are derived for the spin-singlet BB states in 
accordance with the required parity conservation. 

The complete transition matrix elements for charmonium decay into BB states 
(with spin projection /z) can be written as 

M„ = E f« V «( B ) clmI^L^(i-m) (4) 

q.LM 

where the parameter f q represents the amplitude of creating the specific qq pair. 
The space integral reads 

4-e = I Jd 3 k <j> L (k) Y LM {k) <j)* B (B - !*) «^(|fc - B) (5) 

with B standing for the final state baryon momentum and <pB denoting the intrinsic 
spatial wave function of the baryon. We have assumed a simple gaussian spatial 
distribution for both the intermediate qq state and for the baryon intrinsic wave 
function. Only the L = and L = 2 partial waves contribute to the qq state and 
we have used 

<j>o (x) 
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= Tl2 V47T 



7-1/2 



3/2 x 2 

eXp( -8^ } 



Q!7r 



1/2 



3/2 



1 X X 

60^ eXp(_ crf ) 
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where the factors til must satisfy the normalization condition n L — 1. As there 
is no evidence that the strength of qq formation process depends on L we simply 
assume no = 712 = -\/l/2 in our calculation. The harmonic oscillator cigenfunctions 
are used for the baryons in their center-of-mass system, e.g. 



P O) = V4tt , exp( 



1 \ 3 / 2 



12 / 1 x 3/2 



2 ' 



(J) 



for the proton and the TV* (1440) resonance. The momentum x stands for the relative 
qq momenta and for the quark-diquark momenta in Eqs.© and 0, respectively. 
The harmonic oscillator parameters a and f3 characterize the size of the relevant 
interaction and we varied them to fit the experimental data. 

It was already noted in Ref.pP that relativistic description seems to be more 
appropriate for the spatial distribution of the final state baryon clusters. Then the 
integral JSJ includes the jacobian due to Lorentz transformation from the baryon 
CMS to the laboratory (J/tp at rest) system and the internal quark-diquark mo- 
menta are transformed appropriately as well (see 012] for more details). We have 
performed the calculation for both the nonrelativistic as well as the relativistic 
baryon wave functions. 

The space integrals (0 were computed numerically by using the program RI- 
WIAD from CERN Program Library. The decay cross-section for J/tp — > BB was 
then constructed from the amplitudes M„ as 



dT(J/tp BB) 1 I B 



E 1 M » 1 2 ( 8 ) 



dO, 32tt 2 Ml 

with denoting the solid angle of B and M^ standing for the charmonium mass. 



3 Results and discussion 



There are three parameters in our calculation. The harmonic-oscillator param- 
eters a and /3 determine the shape of the angular distribution 

dTU/ib BB) , rj 

V ' dn ! -=N SB (l + a B cos 2 9) , (9) 

i.e. the computed parameter a B ■ The constant Np p is directly related to the experi- 
mental branching ratio of J/ip — > pp and can be used to fix the overall normalization 
of the computed rates. Finally, we define our third parameter as the rate g = f s /f u 
and naturally assume fd = /„. Clearly, the rate g affects the relative numbers 
of strange and nonstrange baryons produced in J/tp decay. As the creation of qq 
pair in charmonium decay is flavor blind one may simply assume g = 1. However, 
the variation of g also accounts effectively for some differences between creation of 
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Table 2. The computed characteristics Fg B /Fp p an d a B of the J/ip — * BB decay rates. 





nonrelativistic case 


relativistic 


case 


experiment 


BB state 


^bb/^pp 


a B 


^bb/^pp 


a B 


^bb/^pp 


a B 


PP 


1.00 


0.64 


1.00 


0.60 


1.00 


0.61(11) 6 


pN* 


0.90 


0.53 


0.76 


0.68 






N*N* 


1.05 


0.18 


0.99 


0.32 






AA 


0.56 


0.50 


0.54 


0.54 


0.61(9)13 


0.62(22)0 


E°E° 


0.41 


0.50 


0.39 


0.56 


0.60(11)0 


0.22(31) 6 



various DD species, i.e. between different amplitudes fj^ D - In the calculation we 
adopted the value g — 0.9 suppressing partly the formation of strange baryons in 
the final state. Our assumptions and the structure of the amplitudes V q (B) allow 
to sum over q in Eq.@. The summed amplitudes J2 q fqVq(B) are gi yen m the last 
column of the Table [5] 

The results of our calculation are presented in the Table for both the nonrela- 
tivistic and the relativistic approaches. We show the relative decay rates T^ B /Tp p , 
where Tg B = Ai:Ng B {l 4- as/3), and the angular distribution coefficients as in 
comparison with the available experimental data. As the reliable data are limited 
only to the pp channel we futher assumed a = (3 and made an one parameter fit 
to the measured value of a p . The results shown in the Table were obtained for 
a = (3 = 0.4 GeV and for a — (3 — 0.22 GeV in the nonrelativistic and the relativis- 
tic cases, respectively. Although the assumption of using the same size parameters 
for both the intermediate qq state and for the space distribution of baryon clusters 
may not be sound the fitted values compare well with those used in other quark 
models 00. 

The results are compatible with available experimental data on pp, AA and 
S°E channels within two standard deviations. However the predicted results for 
the pAf*(1440) and N*N* are quite different from those given by the simple quark 
model 1 which predicts rates for these two channels that are at least two times 
larger than the values presented here. Future experimental results on these two 
channels will be helpful for examining various model predictions and to improve 
our understanding of internal quark structure of these baryons. If the experimental 
N* production rates turn out to be much lower than the quark (and diquark) 
model predictions a large component of 7rA^ could contribute to the N* internal 
structure. On the other hand, too high production rates would indicate some other 
mechanisms playing a role. 
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edge Innovation Project (KJCX2-SW-N02) and the National Natural Science Foundation 
of China. A.C. acknowledges the hospitality of the IHEP in Beijing. 
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Appendix 

In this section we show the quark-diquark wave functions of baryons used in 
the present work. Suppressing color anti-symmetrization factors the spin-up and 
spin-down proton states read 

| p±) = -^-j= |±[ud] u± T V2[uu]od± =F V2[ud]±\u T ± 2[uu]± 1 d T ^ + -^=(ud)u± . 

(10) 

The N* resonance has exactly the same spin-flavour structure and we can write 
for the A and S°: 

| A±) = ^±[ds} u± T [us] d± T V2[ds] ±1 u T ± V2[us] ±1 d T + 

+ (ds)u± + (us)d± - 2(ud)s±} (11) 

| Si) = ^{±2[ud} s±T[us}od±T[ds} u±T2V2[ud} ±1 s T ±V2[us} ±1 d T ± 
± V2[ds] ±1 u T } - - {(us)d ± + {ds)u±} (12) 
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